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Abstract. There are three quarks with masses at or below the characteristic scale of QCD dynamics: 
up, down and strange. However, twisted mass lattice QCD relies on quark doublets. Various options for 
including three quark flavors within the twisted mass approach are explored by studying the kaon masses, 
both analytically (through chiral Lagrangians) and numerically (through lattice simulations). Advantages 
and disadvantages are revealed for each "strange and twisted" option. 
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1 Two flavor twisted mass lattice QCD 

Consider a hypercubic lattice in four-dimensional space- 
time, and a two flavor system: 

1 *■>-($)• (i) 

In the m u = rrid = m q limit, the fermion action for twisted 
mass lattice QCD is 

x 

, with Dirac and Wilson terms, 

D(r,uj) = 7 - V + exp(-tW7fc73i)W(r) , (3) 

(if 

W(r)^(x) = M cr (r)^(x) - — □ ^(z) . (4) 

As usual, the covariant derivative and dAlembertian 
□ contain the gauge field U^x) = exp(iagT b A b Jx)) . 

Two key features of twisted mass lattice QCD are (a) 
the existence of a firm lower bound for the eigenvalues of 
D(r,u>) + m q , as long as m q and u> are non-zero|I] and 
(b) the absence of 0(a) errors in observable quantities 
at maximal twist, uj(renormalized) = ±7r/2[2]. Another 
significant consequence of w ^ is that parity and fla- 
vor symmetries are broken; they only get restored in the 
continuum limit. 

2 Options for including the strange quark 

Given the phenomcnological importance of u, d, s physics, 
how should the strange quark be added to the two flavor 
twisted mass theory of the previous section? 
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One option is to leave the strange quark untwisted. It 
will not be protected from near-zero modes, but this is not 
a practical problem for simulations at the physical strange 
quark mass. There will be 0(a) errors in the theory, but 
one could remove their effects on observables through ad- 
dition of a Sheikholeslami-Wohlert termf^]. Kaon opera- 
tors, containing one twisted and one untwisted fermion, 
would appear as parity mixtures, requiring extra effort 
for the extraction of continuum physics. 

Another option is to twist the strange quark by using 
the charm quark as its twisting partner. To accommodate 
m c ^ m s , the action for this "heavy" doublet will need 
one additional term relative to the "light" doublet action 
used in the previous section for up and down. The quark 
mass terms in the twisted basis are 

C = i>i {mi + ij 5 [iiT 3 ) ipi 

+ iph (mh + ij5l^hT3 + eT a ) ip h , (5) 

where r a is a Pauli matrix in flavor space. Should we 
choose r a — r 3 ^4 or r a ^ r 3 5 ? An advantage of the par- 
allel choice, i.e. r a = T3, is that flavors do not mix, but an 
important disadvantage is that the fermion determinant 
is not real. Therefore this parallel choice corresponds to 
fermions that can be used as valence quarks but not as sea 
quarks. The perpendicular choice, r a ^ T3, leads to a real 
fermion determinant so those quarks can be both sea and 
valence, but the flavors in the action mix, as is evident 
from |J5J|. To elucidate the flavor and isospin structure of 
this choice, we study the kaon masses in this work. 

The notation (c, s)y and (c, s)j_ will respectively de- 
note twisted doublets with the parallel and perpendicular 
choices for r a relative to T3. A degenerate doublet is de- 
noted by (u,d)o. The three scenarios to be explored are 

scenario 1: (u, d)o + s [+c if desired] , 

scenario 2: (u, d)o + (c, s)ii , 

scenario 3: (u, d)o + (c, s)± . 
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3 Kaon mass splittings and lattice artifacts 
3.1 Chiral perturbation theory 

The chiral Lagrangian for two twisted doublets is(H] 

£ xPT = £^ + ... , (6) 

£(2) = Lt^d SD tf) _ l-Tiix^ + E^x) 



Tr(A* E + A) 



£ (4) = -W£Tr[(ti £ + A) 2 } + 



where 

£ = 



,lWZT3/2 





e iuJiT 3 /2 







,iuJ h r 3 /2 



(7) 
(8) 

(9) 



for meson matrix <!>. The quark mass matrix appears as 
X = 2 B ( mi+ ^ , . , V (10) 

and the lattice spacing as A — 2Woffl, using notation fa- 
miliar from |7I5I5| . 

To interpret (u, d)o + (c, s)j_, we diagonalize the quark 
mass matrices, 



(11) 



m h + ifi h r 3 + e hTl = e ^3^/2 F t M y e mc h /2 ) {u) 
where Y 



l 

v/2 



\ Jj 1S a standard 7r/4 rotation, and 
the diagonal matrix M is 



M 











(13) 



The resulting mass terms for kaons and D mesons have 
the form 

(K+ ©° K° ED" 








) 









( 14 ) 

where m, (5 and a contain various parameters from the chi- 
ral Lagrangian. Here, mesons are labelled by their iden- 
tities in the untwisted limit. They are clearly not mass 
eigenstates for nonzero twist, a ^ 0. To determine meson 
eigenstates for a general twist, we diagonalize the matrix 
in i|14|) through a simple change of basis, and obtain the 
eigenstates and eigenvalues of Table ^ The rotation (not 
twist!) angle, < 9 < tt/2, is given by 



tan 29 



(15) 



where a ~ sin lji sin uj h vanishes if either [u, d) or (c, s) 
is not twisted, and 5 ~ vanishes when (c, s) has no 
explicit mass splitting term. 



Table 1. Eigenvalues and eigenstates of the matrix in l|14fl . 
Note that, for general 6, none of these mass eigenstates are 
isospin partners. 



mass eigenstates 



tpu- = K + cos 6 - D u sin ( 
ipd- =K°cos6> + B~sin( 
ip u+ = K+ sin + D° cos ( 
ipd+ =K°sin6l-D _ cos( 



mass eigenvalues 



m u - = tnd- ~ m — \J5 2 + a 2 



m u+ = m d + = m + VS 2 + a 2 



Is there a correspondence between these eigenstates 
and the physical mesons? As a special case, consider a = 
0. This implies 9 = and leads to scenario 1: (u, d)o+s+c. 
Quark flavors do not mix; the physical kaons are a (lighter) 
degenerate isospin doublet, and the physical D mesons are 
a (heavier) degenerate isospin doublet. This scenario is 
also discussed in [TU| . 

As a different special case, consider 5 = 0. This implies 
9 = 7r/4 and leads to (u, d)o + (c, s)o- Quark flavors can 
now be diagonalized, and mixings are thus avoided. In 
fact, this special case then becomes identical to (u, d)o + 
(c,s) II with m c = m s . The mass-degenerate doublets are 
(K + , D~ ) and (K°, D ). The physical kaons only become 
degenerate in the continuum limit. 

In general, when a ^ and 8^0, Table ^ corre- 
sponds to (u, d)o + (c,s)j_. The light quark doublet has 
no flavor mixing, but the heavy doublet does. Interest- 
ingly, none of the mass eigenstates form isospin doublets 
unless isospin is somehow defined to involve (c, s) as well 
as (u,d). Our findings are sketched in Fig. [2 Scenario 1 
has no isospin splittings among kaons, and none among 
D mesons either, but there are 0(a) splittings between 
the kaons and D mesons. In scenario 2 all splittings begin 
at 0(a 2 ), but twist artifacts do cause unphysical 0(a 2 ) 
isospin splittings between charged and neutral kaons, and 
also between charged and neutral D mesons. In scenario 3, 
flavors mix. There are two degenerate pairs of eigenstates, 
but the symmetry relating the degenerate states involves 
both quark doublets. Identification of these states with 
physical particles only becomes clear upon extrapolation 
to the continuum limit. 



3.2 Lattice QCD simulations 

According to the chiral Lagrangian discussion above, the 
splitting m 2 (K° ) - m 2 (K ± ) that arises from twist artifacts 
in the (u, d)o + (c, s)h scenario should be, at leading order, 
(i) independent of quark mass, (ii) linear in a 2 , and (iii) 
vanishing as a 2 — > 0. Quenched numerical simulations for 
this scenario |dTTiT] are consistent with these expectations, 
though the vanishing of the mass splitting as a 2 — > is 
not as clear in the data as one might have hoped. Even at 
the smallest lattice spacing, a w 0.068 fm, the kaon mass 
splitting is found to be sizeable: m(K°) — m(K ± ) ~ 50 
MeV. 

Ref. [H] also contains quenched results for the (it, d)o+s 
scenario, and these confirm the chiral Lagrangian claim of 
m(K°) = m(K ± ), i.e. the kaon masses are not split by 
twist artifacts. 
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Fig. 1. Meson mass eigenstates obtained from three scenarios, 
(a) scenario 1: (u,d)o + s + c, (b) scenario 2: (u,d)o + (c,s)ii, 
(c) scenario 3: (u, d)o + (c, s)±. 



For scenario 3, (u,<i)o + (c, s)j_, conventional isospin 
is only restored in the continuum limit. Therefore meson 
mass eigenstates must be identified by computing the full 
meson correlation matrix and then diagonalizing it numer- 
ically. With quark mass terms 



4 Summary 

Twisted mass lattice QCD removes 0(a) errors and un- 
physical near-zero modes. The up and down quarks fit nat- 
ually into a twisted doublet but the strange quark has no 
natural partner. Three scenarios for including the strange 
quark were considered within chiral perturbation theory, 
and found to have the following features: 

• (u, d) + s + c 

— no unphysical mass splittings within isospin multiplets, 

— the strange quark does not benefit from twisting, 

— kaon operators (i.e. twisted+untwisted) require care, 

• (u,d) + (c, s)|| 

— all quarks benefit from twisting, 

— suitable for valence quarks but not for sea quarks, 

— flavors do not mix and isospin is easily managed, 

— discretization effects appear as kaon mass splittings, 

• (u,d) + (c, s)_l 

— all quarks benefit from twisting, 

— dynamical simulations can use a single action, 

— flavors mix; isospin is not straightforward for a ^ 0. 

These features have also been observed in explicit numer- 
ical simulations, and may be helpful when choosing which 
scenario is optimal for a particular phcnomcnological ap- 
plication of lattice QCD. 
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(u d) ( mi + ij 5 [ii 

^ mi - z7 5 /i; 

(q~i q 2 ) fm h + i^h e h 



e h m h - ij 5 Hh I \12 



(16) 



and operators U i — wy^qt and Odi = d-y^qi, numerical 
results give 



cp = (o ui ol)^c^ = 



a c 
c b 



(17) 



(°*°l)=^« = (^), (18) 



where a, b and c are real-valued. Diagonalization of 
and provides the meson mass eigenstates. Simulation 
results [T21ll3| are consistent with our chiral Lagrangian 
discussion. 
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